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Abstract 


In eikonal and quenched approximation, it is argued that the strong coupling fermionic QCD 
Green’s functions and related amplitudes depart from a sole dependence on the SUc{3) quadratic 
Casimir operator, C 2 /, evaluated over the fundamental gauge group representation. 

Noticed in non-relativistic Quark Models and in a non-perturbative generalization of the 
Schwinger mechanism, an additional dependence on the cubic Casimir operator shows up, in 
contradistinction with perturbation theory and other non-perturbative approaches. However, it 
accounts for the full algebraic content of the rank-2 Lie algebra of SUc{3). Though numerically 
sub-leading effects, cubic Casimir dependences, here and elsewhere, appear to be a signature of the 
non-perturbative fermonic sector of QCD. 
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I. INTRODUCTION 


In some recent articles [l-5|, a property, which bears on the non-perturbative fermionic 
Green’s fnnctions of QCD, has been put forth under the name of effective locality. This 
property can be phrased as follows. For any Quark/Quark (or Anti-Quark) scattering am¬ 
plitude, the full gauge-invariant sum of cubic and quartic gluonic interactions, fermionic 
loops included, results in a local contact-type interaction, and this local interaction is me¬ 
diated by a tensorial held which is antisymmetric both in Lorentz and color indices. This 
is a non-expected result because, ordinarily, integrations of elementary degrees of freedom 
result in highly non-local structures; the ‘effective locality’ denomination, which sounds 
like an oxymoron, accounts for this rather unusual circumstance. It is worth emphasizing 
that effective locality is an exact property of full QCD, and that its derivation entails no 
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approximation 

Then the consequences of effective locality, even when examined ‘at tree level’, should 
exhibit admissible as well as new aspects of the conhned phase of QCD; and so far, it seems 


5 |. 


to be so 

In Ref. , a general form of the non-perturbative QCD fermionic amplitudes is displayed 
as a hnite sum of hnite products of Meijer special functions , in agreement with general expec¬ 
tations jb]. Remarkably enough, within one and the same expression, these amplitudes are 
able to display an explicit link between a partonic content and a hadronic non-perturbative 
component in accord, this time, with the AdS^/QCD light-cone approach of Ref.[^. 

However, the analysis presented in Ref. is carried out at eikonal and quenched approxi¬ 
mations. Soon it will become important to relax these approximations, not only for the sake 
of preserving unitarity, but also in order to explore larger distances: Effective locality, in 
effect, clearly differentiates QCD from a pure Yang Mills situation. In particular, as noticed 
also by lattice approaches, inclusion of quark loops reveals to be essential to the description 
of larger distance non-perturbative physics 

Fortunately, some things can be learnt already at the level of a quenched (and eikonal) 
analysis. In a recent letter js], it was argued that non-perturbative fermionic Green’s func¬ 
tions, and amplitudes thereof, do not only depend on the quadratic Casimir invariant, C 2 , 
but also on the trilinear Casimir invariant, C^. 

This extra dependence obviously complies with the rank-2 character of an S'Dc(3)-color 
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algebra, and at this level of approximations at least, it is a pecnliar ontpnt of effective 


locality. With the exception of Ref-js,!!^ some years ago, it is remarkable that has come 
nnnoticed for so long, whereas its importance was pnt forth in nonrelativistic qnark models 
fn |. Now, the relative smallness of Csj-dependences, as displayed in Section IV, is certainly 
at the origin of this fact. 

Be it as it may, it will certainly matter to disentangle snch a prediction from the approxi¬ 


mated context where it was hrst discovered 


12 | . At present, thongh, the cnrrent paper aims 




at making more complete and accnrate the non-trivial resnlts hrst annonnced in Ref. | 8 | and 
display some nnmerical estimates. 

The paper starts from an expression for the fermonic 4-point fnnction which is the matter 
of a whole article For the sake of traceability. Appendix A offers a snmmary of the main 
steps at the origin of that expression. Besides, the paper is organized as follows. 

For an easier and explicit presentation. Section II introdnces the matter in the case of a 4- 
point fermionic Green’s fnnction, while providing with Appendix D, the necessary extension 
to the general case of 2n-point Green’s fnnctions. In order to alleviate a presentation which 
is already qnite technical, some proofs are deferred to Appendices B, D, E and F. 

Section III deals with the necessary average to be taken over the orthogonal gronp of 
matrices OAr(M), the procednre which reveals the fnll color algebraic strnctnre of 2 ? 7 ,-point 
fermonic Green’s fnnctions and related amplitndes. A proof necessary to Section III is 
prodnced in Appendix G, and a concinsion is proposed in Section V. 


II. NON-PERTURBATIVE FERMIONIC GREEN’S FUNCTIONS 


In Pertnrbation Theory, all of the scattering process calcnlations come ont proportional 
to either Ca = and/or Cp = (A/ — l)/ 2 Ac 5 that is to the qnadratic Gasimir operator 
eigenvalne 6 * 2 ( 7 ^) over the adjoint and fnndamental representations respectively (exceptions 
may be fonnd in the lattice-gange theory approach of Ref. 13|], where higher dimensional rep¬ 
resentation spaces were considered, bnt again, restricted to C' 2 -dependences). The qnadratic 
Gasimir operator’s well known dehnition is C 2 {TZ) = where the Ta{TZ) denote 

the SUc{3) Lie algebra generators in a given representation TZ. 

Not only perturbative calculations, but also non-perturbative QGD models, such as 


the MIT bag model [l^, the Stochastic Vacuum Model (SVM) |l5|, and lattice ap- 
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proaches 161, comply with these overall C 2 {TZ) dependences, though, sometimes, it 
is worth noticing, in quite different ways. It could be that, indeed, perturbative inputs 
be imported as surreptitious elements, from perturbation theory to these non-perturbative 
attempts fl. 

The property of effective locality surfaces at the level of non-perturbative 2n-point 
fermionic Green’s functions, as an exact, non-approximate property of QCD [ 2 ]. In a strong 
coupling regime of g » 1, evaluating a 2 to 2 Quark(Anti-Quark)/Quark(Anti-Quark) 
scattering amplitude with the help of (Analytically continued j^) Random Matrix Theory 
17|, one finds a result proportional to |^, 
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where the eikonal and quenching approximations have been used. QCD is here simplified 
to the case of a single quark species of mass m, and E{= Ei = E 2 ) is each of the two 
colliding quarks energy in the center of mass system, pi = {E,0,0,p), p 2 = {E,0,0, —p), 
■s = {Pi +^ 2 )^- In Appendix D (and thanks also to Appendix F), it is shown how this 4-point 
expression generalizes to the case of 2?T,-point fermionic Green’s functions. 

Note that expression ([1]) is the matter of a whole paper { 4 ]. For the sake of providing 
it with enough background, though, the main steps of its derivation are summarized in 
Appendix A. 

The variables entering ([T]) are the following. 

- In ([T]), O = 0{..,pj,..) stands for an orthogonal N x N matrix specified by the N{N — 
l)/2 parmeters ms, with N = D x (A(? — 1), D, the number of space-time dimensions; that 
is A^ = 32, 24,|25(]. The distribution f{..,pj,..) defines the Haar measure of integration over 
the orthogonal group 07 v(R), and the constant G, its normalization (he., C is the inverse 
of the Oat(M)- volume). 

- In ([T]), the N X N orthogonal matrices O act upon an A^-vector of matrices T = 
(1,1,1,1) 0 T = (T, T, T, T). That is, T is made out of D = 4 copies of the full set T 
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of SUc{3) generators, taken in the fundamental representation: T = with 

ta = Aa/2, the standard Gell-Mann matrices 18|. However technical, this step, related to 
flA7ll . allows one to perform a series of exact integrations in a systematic way. 


- The third line of ([T]) displays a Meijer special function, How these Meijer functions 
come into play is displayed by flAlOjl . Appendix A. The Meijer function depends on an array 
of parameters, one of them involving an integer g*, with 1 < i < iV and 0 < < N{N — 1)/2. 

In Appendix A, Equation (]A9p . it is shown that this power of g^ (see also (lAinp . at p = 
gj —1/2) comes from the expansion of a Vandermonde determinant into a sum of monomials. 
It is this sum which is explicitly referred to in the hrst line of ([T]). 

In one and the same argument, it is worth pointing out that the GI 4 Meijer function’s 
argument mixes up partonic variables (m, s) with the non-perturbative function. 




( 2 ) 


where b = \b\ is the impact parameter of the 2 scattering quarks in their center of mass 
system, and p, the mass scale necessarily introduced by the property of effective locality jj], 
as stated in Appendix A, after Equation flAbl) . This is interesting because, as quoted in the 
Introduction, that point turns out to be a prediction of the AdS^/QCD approach to the 


the order of one tenth j^. While this parameter is crucial to a description of conhnement 
in terms of a linear potential between quarks [l-4|, it has no relevance at all to the current 
considerations. Likewise, the non-perturbative function of (Ej) to be discussed in details else¬ 
where, is partly phenomenological; this point, however, has only a very marginal incidence 
on the results of the present paper. 


non-perturbative regimejDf QCD [^. In ([2]), stands for a small deformation parameter, on 


Two remarks are in order. 

(i) Deriving ([1]), the absolute values of flA 8 j) . relevant to the QCD case of k = 1, have 
been dropped. That this simplihcation can be made without compromising the overall color 
algebraic structure of Green’s functions, is shown in Appendix B. 

(ii) The color algebraic part of (II]) appears in its 2 nd line in terms of the 2 sine functions 
with their arguments being the W x A/ matrices [ay{OT)i], J = {1, 2 }. For each monomial 
of the Vandermonde determinant flA9jl . the structure is that of a hnite product of terms 
attached, each, to a given value of the index i, ranging between 1 and N, and no summation 


5 













understood on i. At this order of approximations at least, this amounts to the statement that 
fermonic Green’s functions and related amplitudes split into a hnite sum of hnite products 
of Meijer special functions. 

As will be further commented, this simple structure is far from being a trivial point, as it 
could be connected to the deep meaning of effective locality. In Appendix F it is shown that 
Baker-Campbell-Hausdorf formulae would dehnitely object to it. However, this structure 
emerges as a consequence of the On(M) average to be taken in ([1]). 


Now, as shown in Appendix E, it is possible to rely on a series of identities and exact 
textbook integration formulae, so as to take ([T]) to the form, 
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Ojv(IK) 


where the GqH Meijer function is explicitly dehned in flAlOp . In ([3]), the large brackets are 
here to denote the orthogonal group 07 v(K)-average specihed in the 2nd line of ([1]). In 
random matrix calculations, it often happens that the eigenvalue spectrum only matters, 
and that the (9Ar(M) averages just factor out and disappear in the normalization. This is 


not the case here: As pointed out in Refs. [4, l8| in effect, the additional but unavoidable 
complexity coming from the (9Ar(M)-average is essential not only to prevent a trivial result 
from occuring, but also to display the full algebraic content of the rank-2 SUciff) color 
algebra. 

The matrix-valued argument of the Meijer’s function of (|3]) is. 


z, = X[{OT)f\\ 


A = 


^ ( 9^{^) 
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(4) 


Y y/32iNc 

As argued in Appendix E, even at a sufficiently large value of the coupling constant g, the 
matrix elements of Zi are much smaller than unity, each. It is therefore possible to rely on the 
standard Meijer function expansion, flE3jl . flE4j) . and write for a given monomial appearing 
in the sum (|3]), 

N 3 
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where the bh stand for the parameters of the G-Meijer function of ([3]), that is bh = 
{bi, & 2 , ^ 3 } = {1/2, (3 + 2qi)/A, 1}. Out of Appendix E, numbers A/ are readily identihed to 
be given by, 


A—r(-)r( ^ 




and likewise. 


Oil — 


), A = r(- 


Oio — 


^)r(- 


^ ), 4 = r(-i)r(A2). (6) 


16 


2qi-3' 


(2g, + 3)(2g, + 5)’ 


Oi3 — 


2qi + 3 


(7) 


Note that, with no incidence at all for the sequel, (E]) brings a correction to Eq.(8) in Ref. js]. 


III. 07v(m)-integration and casimir operator dependences 

Thanks to effective locality, the original (and inhnite dimensional) x“j,-functional sum¬ 
mations of flA2p . 

8 3 

f dki = n n n / <iKM] («) 

wGA 4 a=l 0=/i<iy P 

can be translated into an analytically continued Random Matrix integration ^ (in (|8]), 
A4 is the four dimensional spacetime manifold). The latter splits into an integration on 
the spectrum of matrices M dehned in flA7p . and an integration over the orthogonal group 
As displayed in flAlOp . the former yields the Meijer functions, whereas the full color 
algebraic dependences of fermionic Green’s functions come about as the OAr(M)-integration 
is carried through. 

The matrix-valued argument Zt is on the order of [(C>T)A, and thus, the orders z^, yA, 
and Zi^pZi contributions of ([5]), are of even orders (C>T)°, {OT)1, {OT)j and iOT)^ respec¬ 
tively. They are leading contributions in view of the smallness of A (Appendix E), whereas 
odd powers of (C>T)iS, such as 2qi + 1, vanish trivially under OAr(M)-averaging. Random 
orthogonal matrices, O, can be generated in different ways, distributed according to the 
Haar measure over the orthogonal group (9Ar(M) 19|]. An orthogonal matrix is conveniently 
decomposed into a product of N{N — l)/2 rotators, plus reflections. 


O — (R 12 R 13 ■ • • Rin) (R 23 -R 24 • • • R 2 n) . {Rn- 1 ,n) De , 


(9) 


where the matrix of reflections is diagonal by definition and reads, = diag{ei, 82 , ■■■ ,Sn), 
with Ei = ±1, Vz = 1,..., A. A random orthogonal matrix requires that to either value 
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Si = ±1 an equal probability of P{ei = ± 1 ) = 1/2 be associated. A rotator Rij{Qij), itself 
an X A^-orthogonal matrix, acts as a rotation in the (i — j)-2-plane solely, and is thus 
characterized by an angle 0 jj, while being restricted to the identity operator 12 x 2 on any of 
the other (Z—m)-2-planes with either l,m ^ i or l,m ^ j. The 0,^ are independent random 
variables with a joint probability distribution proportional to 19|, 


N 

n 

i=2 


N 

COS-^~^ 01j 


N 


COS' 


i-3 


02j ■ ■ ■ n ^ ’ 


1=3 


j=N 


whereas the probability density of an angle Qij is a beta distribution, f, with 

cos0ij = yTij, that is, fi{xij] |, |) = x^-"^ ^(1 — |). As meant in the second 

line of ([T]) , these probability densities allow one to calculate averages over orthogonal matrices 
in a dehnite quantitative way, though somewhat probability density dependent. 


For our purpose however, this unwanted probability density dependence will not affect 
our derivations as the full explicit form of the Haar measure is not required, but only the 
Tb matrix properties, and the left- and right- invariances of the Haar measure on Ojv(lR) 
JJI]. Denoting by .., Qim ■,■■■) the matrix elements of ([9]) as the reflection matrix 

Dg is omitted, one obtains. 


= {ep\Q) = S^a^^P)a^\Q))JP , ( 10 ) 

where, in the last equality, the average over the product of reflections ejSk has been taken. 
That is, < pZi > is given by 

N N /-r 

j,k=l 1=1 

where 6 * 2 / stands for the quadratic Casimir operator eigenvalue on the fundamental repre¬ 
sentation, 6 * 2 / = Cf = 4/3, and where Isxs is the identity matrix of format 3x3. The 
last equality of flTTD . independent of the index i, is a consequence of the left- and right- 
invariances of the Haar measure on OAr(M), as it is shown in Appendix C . 


^a|(... Qim ■ ■ ■ ))qT^ - DC2f Isxs , (11) 



In the same way, for the sub-leading piece of ([5]), one obtains, for alH = 1,..., A^, 




in which one notices the occurrence of the cubic Casiniir operator eigenvalue in the funda¬ 
mental representation of S'f/c(3), 

N^-l 

= C^f 1^x3 ■ (13) 

a,6,c=l 


The fully symmetric constants dabc are dehned in the standard fashion, that is, in the case of 
interest, for = 3, {4,4} = dabctc + l^ab- Far less popular than C 2 , the trilinear Casimir 
operator eigenvalue over a representation space specihed by the Young Tableaux parameters 
(p, q) is given by 


C3{p,q) = Y^(p-g)(2p + g + 3)(2g + p + 3), 


(14) 


and Csf = 43 ( 1 , 0 ) = 10/9 in the SUc{3) fundamental representation [20|, whereas it is 


C^a = 43 ( 1 ,1) = 0, in the adjoint representation, (1,1), another salient feature of distinction 
between QCD and the pure Yang-Mills case. 


At next to sub-leading order, Zi^pZi^ corresponding to the Oat(M)- averaged value of 
(4*l(p)7))®, calculations become more intricate. One hnds, 

{ (2(£>C2,)2 + (DC2,)(DC,f) + l(r>Cv)) 13X3 

+ ^2 dkjmdkhi {TjTmThTi + 2TjThTiTjn)\ ■ (15) 

k,j,l,h,m ^ 


While independent of i, the two last terms look somewhat puzzling as they seem to com¬ 
promise the general structure of these dependences, in terms of SUc{3) algebraic invariants. 
To proceed, one may rely on the standard textbook values of the dabc coefficients [18|, and 
work out the following two identities. 


dkn-O, J2 dk'jmdkjm- ^dk^k 

k,j=l j,m=l 


( 16 ) 


of which, the second one can also be found in Ref. 2l|. It is then easy to prove that (1T5|1 
indeed reduces to. 


< Zi^/zl > — (-^)^ ^ [2 + (g)^](-D42/)^ -h {DC2f){DC5f) + 3(1143/)^ 13 x 3 • (17) 

An analysis of the higher order terms of ([3]) is likely to become much more complicated and 
not very insightful either, as it is now made clear enough that 43/-dependences will enter 
any 2n-point non-perturbative fermionic Green’s function and related amplitude, and this 
is the point of the present paper. 
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IV. SOME NUMERICAL INSIGHTS 


Concerning numerical orders of magnitude, the result of flT^ shows that at sub-leading 
order (\/A/iV)^, the trilinear Casimir operator C^f enhances the pure Cb/ contribution an 
amount of 15.6%, to be compared to the 15% at most, advocated in [lOt], whereas at sub- 
sub-leading order (y/X/N)^, C 2 / and Csf contributions to ffT7|) are identical to within 0.2%. 


Now, even at a very large absolute value of the A-parameter of (jl]), (see Appendix E, 
Eq. flE2p ). the contribution brought about by the Cs-dependence enhances the pure C 2 / 
contribution a small relative amount of, 

%/A 


Da 


3/ 


~ 0.01%, 


(18) 


N DC,f + (^) {DCyY 

where y/X/N ~ 0.012%, a value obtained at a strong coupling of = 15, 

/i^/sexp — ~ 1, m = 5MeV and s = lOOMeU. Likewise, as compared to a pure 

linear dependence in C 2 / alone, one gets a less than 0.1% relative deviation, 

%/A iDC2fy + DC3f 


~ 0.075% 


(19) 


N DC2f 

Though it is here dealt with dynamical rather than static quarks, so small departures 
rom a pure C 2 /-dependence comply with the roughly measured linear C 2 /-dependence of 
3 2nd paper), as well as it supports the ex per imental 5% of maximal deviation from the 
so-called C 2 f-scaling hypothesis advocated in [l3| (1st paper). 

The point however, is to estimate how the elementary relative deviations ffTSj) or flT^ 
translate at the level of a whole monomial, and, further, to the full sum of them. This is 
quite long a numerical affair, two hints of which may be given here. 

- Neglecting the sub-sub-leading correction flTTl) to the leading and sub-leading ones, 
(fTOj) and (ITT]) respectively, and considering a monomial like ^^^ 2 ^ ■ ■ ■ associated to a 
distribution of the powers N — 1, N — 2, ■ ■ ■ , 1, 0 for the qiS, one gets for (j5]) the expansion. 
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where ([6]). ([7]) as well as the relation r( — i)/r(i) = —2 with r(|) = ^/^T have been used. 
Introducing the shorthand notations, 

-2r(^) 


. '’•- 25 ,- 3 ’ 


( 21 ) 


equation ([T9]) can be reported to the level of a whole monomial, giving, as a relative deviation 
5c, to a pure C' 2 /-hnear Casimir behavior. 


5c = 


, nri ^ Ef + E?" a* El ^ 


+ DCy) 


+ o{{- 


( 22 ) 


N'-Cy' 

For the monomial considered, involving the powers iV — 1, iV — 2, • • • , 1, 0, a long calculation 
yields a relative departure of. 


5c(iV - 1, iV - 2, • • • , 1,0) ~ -1.0% , 


(23) 


which isn’t that small a number if we keep in mind that there is an enormous number of 
monomials. However, a large number of monomials involve the same values for the g^s, 
attributed to different eigenvalues ^iS to be integrated upon. Alternate in + and — signs, 
these contributions to the full sum of monomials cancel out. This is made explicit within 
the following example, worked out in the easiest case of = 4. 

- At = 4, one has 2® = 64 monomials, that reduce to 48 due to a trivial cancellation. 
Each monomial is affected with a plus or a minus sign and these signs 

come about in an equal number of occurrences. Once integrated over the full spectrum of 
eigenvalues, the and the auxiliary field variables, (lAlOp and (lAllll . the contribution of 
a given monomial can be represented by a ‘word’, 

±eer---e»r ^ ±(gig2---gw), ( 24 ) 

with the sum of the q[s adding up to N{N — l)/2 = 6, in this case. In this way, the sum of 
all 48 fully integrated monomials admits the following symbolic representation, 

+ (3210) - (3201) + (3102) - (3120) + (3021) - (3012) - (2211) + (2202) 

-(2103) + (2121) - (2022) + (2013) - (2220) + (2202) - (2112) + (2130) 

-(2031) + (2022) + (1221) - (1212) + (1113) - (1131) + (1032) - (1023) 

-(2310) + (2301) - (2202) + (2220) - (2121) + (2112) + (1311) - (1302) 

+ (1203) - (1221) + (1122) - (1113) + (1230) - (1221) + (1122) - (1140) 

+ (1041) - (1032) - (0231) + (0222) - (0123) + (0141) - (0042) + (0033). (25) 
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Averaged over 07 v(®), the sum of contributions (1^ yields ([3]) taken at At = 4, and it can 
be further expanded along fl20|) in view of the smallnessss of \/X/N ~ 0.1% at At = 4. 
Now, clearly, the order along which the numbers qiS appear in a ‘word’ is irrelevant to the 
result: Contributions such as (3210), (1203), (2031), for example, are one and the same 
contribution. Accordingly, a lot of cancellations take place and fl25ll eventually reduces to, 

2 X (0222) - 2 X (0123) - (1212) + (0141) - (0042) + (0033). (26) 


As compared to a purely linear behavior of the fermionic 4-point function in 6 * 2 / alone, 
the extra contributions in and C^f are therefore in a ratio of, 

1798,211444 \/A DCsf + {DC2ff 


= -7. 


(27) 


147,543465 4 DC's/ 

Such a result doesn’t comply exactly with the allowed experimental boundaries of 5%. 
However, several points must be taken into account: First, it holds at At = 4, away from 
the realistic value of At = 32. In a second place, it is obtained within the quenching and 
eikonal approximation. Thirdly, it is obtained on the basis of /r^/s'exp —(/i 6 )^“^ ~ 1, a 
phenomenological input which is here deliberatly overestimated, given that /r the effective 
locality scale is not known at present. In effect, non-perturbative physics should take place 
at a distance b > 1/Aqcd, so that the exponential of exp — (p,&)^“^ could reduce the result in 
a significant way if the effective locality scale is larger than Aqcd (for example, such would 
be the case, would the effective locality scale u be close to some chiral symmetry breaking 
estimates, a possibility under consideration |23|). 

Clearly, a way has to be found to cope with the realistic value of A^ = 32; work in this 
direction is in progress where, to begin with, the value of A^ = 32 has already been reduced 
to the smaller value of A^ = 16. 


V. CONCLUSION 

The newly discovered property of effective locality is here explored in one of its con¬ 
sequences, namely, the color algebraic structure of the fermionic strong coupling Green’s 
functions and amplitudes thereof. At variance with a whole series of non-perturbative ap¬ 
proaches, the latter are found to display dependences on the full color algebraic content of 
the S'7/c(3)-algebra, which is of rank 2. That is, not only G 2 /, the quadratic Casimir operator 
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in the fundamental representation of SUc{3) is involved, but also the trilinear one C^f. In 


the absence of a superselection rule that won 
in a way or another, such a result is natural 

At the exception of a few previous results 


d prevent C^f to show up, one may think that 
16, 


10 |, where similar quadratic and trilinear 


color-invariants were put forward as relevant to the description of non-perturbative calcula¬ 
tions, this is quite unexpected a result and this is why it matters to present some details of 
its derivation. 

This is the purpose of the current article where Appendices are used to keep a main text 
from being overwhelmed with too many technical intricacies. In Appendix A, for example, 
one summarizes the context in which the starting expression ([T]) of the current analysis 
is derived, whereas in Appendix F, one can realize that the welcome structure of strong 
coupling fermionic Green’s functions, in terms of hnite sums of hnite products of Meijer 
functions is a non-trivial result: It is induced by OAr(M)-averaging. Over exponentials of 
non commuting SUc{3) generators, OAr(M)-averaging operates in exactly the same formal 
way as Time-ordering does on exponentials of non-commuting field operators. We think 
that these structures shed some new lights on the genuine nature of effective locality, as 
should be explained elsewhere. 

As it comes to numbers, our preliminary results seem to fall into the experimentally 
allowed range of possible deviations from a sole quadratic Casimir dependence of non- 
perturbative fermonic Green’s functions. Cubic Casimir contributions appear only as sub¬ 
leading effects. On a purely theoretical point of view though, a comparison of Equations 
cm and flT^ shows that a pure linear C' 2 /-dependence is but an approximation, however 
excellent. These results are derived within two major approximations that are eikonal and 
quenching. While the former one is appropriate to describe high energy scattering processes, 
the latter doesn’t preserve unitarity. 

First insights [l2|, obtained by relaxing the quenching approximation display the kind of 
compensation mechanisms that are peculiar to the unitarity constraint. But, even though 
modihed, G 3 /-dependences appear to persist under unitarity restoration. Of course, a more 
quantitative determination of these dependences is in order, and will be addressed in a 
forthcoming publication [l2|. Though subleading, it is quite likely that G 3 dependences be 
the hallmarks of genuine non-perturbative calculations. 
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Appendix A: Context of Equation ([!]) 

Equation ([T]) is non-trivial and is the result of lengthy developments, whose guiding lines 
will be recalled briefly. Relying on standard functional methods, idenBties and notations, 
the QCD generating functional can be brought to the explicit form of 


SqcdIj.i;,!?] f d[x]eUx> fx-r+i I 




where 2)^^ = — | ^ ■ Dp^ ■ ^ is the functional linkage operator, and Dp^ the gluonic 

Feynman propagator in a given covariant gauge with parameter (. The form (lAlIl may 
look a bit unusual as, though equivalent, functional integration is much more customary 
than functional differentiation. Deriving (lAip . y;“j,(x)-Halpern helds have been introduced 
to linearize the original held-strength part of the QCD lagrangian as in 
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(A2) 


The approximation of quenching amounts to take to 0 the fermonic determinant functional 
L[A] of flAip . and redehne the normalization accordingly. Fermonic 2n-point Green’s func¬ 
tions are obtained in the usual way, by differentiating with respect to the fermonic sources 
7 ], fj, and then by cancelling the sources, r] = fj = j = 0. In flAll) . this operation brings down 
functionals G-p{xi,yi\A) which admit exact Schwinger-Fradkin representations, [l-^. One 
has for instance of a mixed representation (space-time and momentum). 


{p\GF[A]\y) = e -i i ds e 


—ism.'" „ —iTrln(2/i) 


• e 2 


(A3) 


X / d[u]{m — zy • [p — gAiy — m(s))]} • e4 


i £ d-s' b'(«')]^ . gipMA 


X 


fo ds'cr-F{y-u{s')) . ^-igf‘ds'u'{s')-A{y-u{s'))'^ 


where s is the Schwinger proper time, u{s), the Fradkin 4-vector held and h the function 
h{si,S 2 ) = /Q*ds'0(si — s')0(s2 — s'). The interesting point is in the last line where only 
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quadratic dependences on the gluonic field at most, appear, but in an s-ordered 

exponential, as indicated by the prescription +. 

The eikonal approximation, then, takes flASp to a simpler form with, to wit, such an 
s-ordered exponential as. 


exp [igpf^ ds' A^^{y - s'p) T“ 


where the T“s stand for the S't/(3)-Lie algebra generators in a suitable representation. By 
introducing 2 subsidiary functional integrations, on a“(s) and ff“(s) fields, say, it is possible 
to take the y4“(a;)-field dependences out of the s-ordered exponential, 

j = J\f J d[a] j d[ff] j 

(A4)' 

and to complete the quadrature in an exact way (i.e., to get the full result of the linkage 
operations that are involved in flAip ). Now, to be guaranteed to deal with the proper 
representation flA3l) or its eikonal approximated form flA4ll . it is essential to perform the 


integrations on subsidiary Q!“(s) and f2“(s) variables in an exact way: This is one of the 
most favourable circumstances which, a few steps later, are met in these calculations jd]. 
Then, the A“(a;)-field quadrature yields the expression, 

where, in the case of a fermonic 4-point Green’s function, 

s“ = -a"xt:„ + 9K„ + Ky. (A6) 

In flA6p . the 7?.“^, i = 1, 2, stand for the eikonal representation involving (IA4p . as the A“(a;)- 
field quadrature has been performed. 


In the_pure Yang Mills case, a form such as 


and 


has been put forth in Ref. 26|. 


In Refs.Q and [8|, as well as in the current paper, the strong coupling limit of >> 1 is 
introduced in a somewhat ‘academic way’, so as to get rid of the V • y [/ • ■ V • y-term 

of (1A5D . that so far, we have not been able to treat on the same footing as the other terms. 


For the reason given just after Eq. llTTll . though, one may argue that this pure gluonic term 
does not impinge on our results. 

Note that the A“(x)-quadrature given above displays the property of effective locality, 
as the generated interactions of the are local, mediated by the structure (/ • x(^)) 


-1 
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taken at the same space-time point z [l|]. The result of the quadrature (besides a residual 
dependence on Dp^ to be absorbed into the normalization constant), does no longer depend 
on the initial covariant gauge-field propagator: The same independence is obtained in a 
similar way with any other bare gauge-held propagator, in any gauge Now, it is of 
utmost importance to emphasize that a series of functional differential identities allows 
one to prove that the same property of effective locality is satished also by the full non- 
approximated QCD theory 2|. This property comes along with a mass scale p, ([2]), (the 
consequence of a theorem whose ultimate identity is still under investigation. 

In the 4-point function case which illustrates the purpose of the current article, and at this 
level of approximation, the inteCTations that are left are those on proper-times Si and S 2 
(made trivial by the theorem of ^), on the subsidiary held variables a“(s) and hl“(s), and 
on the x“j,-helds. 

Concerning the latter, it is convenient to dehne iM, the product, 

W2-l 

® T“ = iM , Mij = Mji e M, 1 < i, j < N , tr M = 0. (A7) 


a=l 


where N = D x — 1) [24, 
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random matrices with measure 
M-i 

T“) 


, so that one has now to deal with an integration over 


nl, 




fJ.l' 


a=l 


= idM = idMii dMi2 • • • dMjvAr 


= i 


d{Mnr-- 

, Mnn) 


■ ■ ,PN{N-1)/2) 


d^l ■ • • d^jv dpi • • • dpAr(Ar_i)/2 


N 


= I 


* n n 1^* “ •• ^pn{n-i)I 2 f{p), (A 8 ) 

i=l i<j 

where the ^iS are the eigenvalues of M, while the pjs complete the parametrization of a 
given matrix M (that is, the original integration on matrices M splits into an integration 
on its spectrum of eigenvalues and on the orthogonal group Oat (M)). In flA8j) . one has a 
Vandermonde determinant of (at k = 1, see Appendix B), 


n 16-01 = 1 E (A9) 

1<2<j<A^ monomials 

The sum in the right hand side of flA9D comprises terms. Each term, a monomial, is 

characterized (not in a unique way) by a given distribution of g^-powers whose sum satishes 
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the constraint of an eqnal global degree of N{N — l)/2. Up to coefficients that are g*- 
dependent, o. o. the monomials share the same algebraic color strnctnre, which is that 
of the net Green’s fnnction and related amplitnde. It is by integrating on a given eigenvalne 


that Meijer fnnctions do appear according to |27|, 

1 


_ _ c 2 _^ 

d^i ef e «» = 


2^/^ 


^03 


p+l 1 


,0 


(AlO) 


In flAlOp . eqnality holds only for p > 0 and 6 > 0, whereas the right hand side is analytic 
both in the argnment and the parameters of the Meijer fnnction: This is the second most 
favorable circnmstance of these calcnlations. 


A third most favorable circumstance comes about as one deals with integrations on the 
subsidiary field variables q:“(s) and G“(s) that have to be carried out in an exact way, as 
quoted above. Defining V/ = OVi, where O is an A^ x A^ orthogonal matrix, and where 
Vi = Pi® G“(0) with the index i = 1,2 labelling the 2 scattering quarks, then the Meijer 
function argument in flAlOp . 6^/4, turns out to be proportional to [(U/)*(V^')*]^. 


Now, in view of (|A4]) or expression and dropping the '-notations for short, these 

dependences on the Vs, must be folded into the following two integrations, 

■ 1 


dU/ e' 


E-p 


■ “2 G 
dU2' e"'W 


20F 


/^30 / “i 
^03 ' ~ 


+ 1 1 p. 

4 ’2’ 


(All) 


28 


29| . Once all pieces are glued togeher. 


whose exact results are provided by the Tables 
the result is given by ([T]) in the main text, the expression the current article is starting from. 


Appendix B: Absolute values in the Vandermonde determinant 


Passing from the functional integrations over the Halpern fields to a summation 
over an algebra of random matrices, one encounters the Jacobian flA8l) . % ll7|, carrying the 
Vandermonde determinant 


n (Bl) 

Now, contrarily to hermitian or symplectic matrices for which the parameter values are 
K = 2 and k = 4 respectively, for a real symmetric matrix, we have k = 1 , and analyticity 
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seems to be lost. But that is not the case. Since matrices M are real traceless symmetric, 
their spectra can be written in the following way, 


Sp M = { (6, ^N-i+i = -6)} , = 1, 2,..., iV/2 . (B2) 


That is, under the form of N/2 pairs of equal and opposite eigenvalues, the eigenvalue 
zero being degenerate, with a mutiplicity greater or equal to two. As a consequence, the 
Vandermonde determinant flBlD boils down to a form, 

N N/2 N/2 


l<i<j i=l 

which, even at k = 1, is obviously analytic. Unfortunately, however, the symmetry prop¬ 


erties of the spectra 


are 


oroken by interaction terms in such a way that the exact 


integration formulae of Ref. 29| can no longer be used. In order to show that the abso¬ 


lute value prescription of (IBip can be dropped without prejudice for our concern, one must 
accordingly construct another argument, more in the line of a ‘physicist proof’. 


Relying on a standard representation of the Heaviside step-function such as 


'• + 0O 


0(a;) = 7 ^ / dz ^ 


2i7r J_^ z — le 


(B4) 


one can write the Vandermonde determinant flBljl i as. 


TT f ^ 


(- 1 )' 


hi pfc+1)! n. 


"+ 0 O „ 2 fc+l 

dz^— 


z — le 


n 


2A:-|-2 


(B5) 


1<2<J<A^ 


Now, in the main text, ([T]) is obtained by integration over the the whole spectrum of fi- 
eigenvalues, with, instead of flBip . the sole term ni<i<j< 7 v subsequently expanded 

into as a sum of monomials. 


n (&-&)= V ± .. . er ■ (B6) 

1<2<j'<A^ monomials 

Bringing flBbp to any power of 2k -|- 2, one gets, 

1<2<j<A^ monomials 

where for every value of fc, the coefficients Cq^{k}...qM{k} as well as the powers qi{k} can be 
written in closed, explicit form. Now, clearly, the coefficients Cq^{k}...qN{k} are irrelevant to 
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the point under consideration, and for the powers qi{k}, the only property of concern is 


29| . that is, V/c, one must have qi{k} > qi, which is obviously 


that it preserves integrability 
satished {qi{k},qi G N). 

With this representation of the absolute value prescription, o should therefore be re¬ 
placed by. 


E 


■1)' 


^H-cxd 


.2/c+l 




d;2- 


z — le 


E c, 

monomials 


qi{k}---qN{k} 


X C dpi .. dpN{N-l)/2 f{pi, . . . ,PNiN-l)/2) 


^+oo 


n 

l<i<N 

sm[a\{OT)i\ I 


da^ 


[1 

i sm[ai{OT)i 


^ + 00 


dal 


a-i 


ah 


xGll 


/ a]al\‘^ s(s — 4771“^) 


2 m^ 


3-2qi{k} 1 1 

4 ) 2’ 


1 

2 ’ 


iff 
2 ’ ’ 


(B8) 


That is, even though a k- series has been introduced, as well as a subsidiary ^-integration, 
the expression fIBSp already displays that all terms share the same algebraic structure which 
is that of the second line of both ([T]) and fIBSIl . This structure can accordingly be analyzed 
on the simpler form flB6p . In full rigor, this conclusion would require that integrations on 
a\ and a^ be carried out, so as to control a possible change of dependences generated by 
the shift qi —)■ qi{k}: That it doesn’t happen to be so is a by-product of Appendix E, 
Eq. flElIl : With qi{k} > q^, replacing g*, the numbers and Oih of ([6]) and ([7]) undergo the 
replacements of A^{qi) —> and Oih{qi) —> Oih{qi{k}), leaving unaffected the (at 

least leading order) Casimir operator dependences resulting from the OAr(M)-average of ([5]). 


Another proof will be given elsewhere 
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Appendix C: O^r (l^)-averages 

As a generic example, one can consider < >e, the OAr(M)-averaged value of 

the product a^j(0)a^^(0), where all 3 indices i,j,k are hxed, and where 0 stands for 
the collective dependences on the angles 0/m, as explicited in ([9]). Let i?/fe(±7r/2) be the 
rotator acting in the I — k plane solely, taking the unit vector > to the unit vector | / >, 
through an appropriate rotation of angle 0/^ = ±7r/2; the restriction of -R/fc(±7r/2) to all 
of the other possible 2-planes is just the unit operator l 2 - Over the whole space spanned 
by the orthonormal basis {|i > ,i = 1,..., |A >}, the standard resolution of unity holds, 

= EJIi K >< ^1- 
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Thus, recalling that Otj = ejatj, with ej = ±1, Vj = 1,..., A^, one has 

< alal>=C f (Cl) 

J O (M) 

where is a standard and shorthand notation for the Haar measure on 07 v(R). Within 
customary Dirac’s notations, fICip can be written as, 

< alal >=C f (< z\0\j >y (< t\0\k >r. (C2) 

Integration on OAr(M) can be right-shifted by multiplying any matrix O to the right, with 
another orthogonal matrix O'] that is, O —?■ OO', 


< 44 >=c [ fiHiOO') (< z\00'\j >f (< z\00'\k >f . (C3) 

By choosing O' = Rik{i:n/2), one obtains 


< 


2 2 


>= C 


„ Af N 

/ iih{ 0) {Y, < AO\h X h\Rik\3 >f CY < X 

dOiv(R) ;j/=l 


\k>) 


= C [ f^HiO)OlOl 

=< 44 > (*^4 


where the right-invariance of the Haar measure is used in the first equality, the definition of 
the particular orthogonal matrix Rik is used in the second one, and fICip in the last equality. 
In the same way, relying on the left-invariance of the Haar measure, the distinguished left side 
index i can be shifted to any other of its possible values in the set{l, 2,..., A^}, establishing 


the uniform probability distribution of these products 


19|. 


From this, and keeping in mind that the aij can be looked upon as the components of 
orthonormalized A^-vectors, it is easy to show (up to some overall multiplicative factors of vr 
which are not relevant to our considerations) that one has < ajj >= l/N, < >= 1/A^^, 

< 444 >= l/^^ etc.. 


Appendix D: Generalization to 2n-point fermionic Green’s functions 

In the case of 2n-point fermionic Green’s functions, one has n-fermionic propagators, that 
is n Fradkin field variables Ui{si), and likewise, n-subsidiary fields ai{si) and Dj(sj). 
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Ignoring renormalization effects for the time being, one has accordingly C‘^ = n{n — l)/2 
2 by 2- interaction terms jd], 

-g / / dsA:nk'^(sfc) / dsiW^{si) 

X -yk + Uk(sfc)) -yi + ui(sz)), 

where 1 < k < I < n. As demonstrated in Ref. jj], each of these 2 by 2-interactions select a 
given point of Minkowski spacetime, Wki, say, where that effective and local k — I interaction 
is taking place, mediated by the strnctnre (/ • Snch Green’s fnnctions therefore 

are proportional to. 


JJ f d^^" f ( 




l<k<l<n ' 
l<a<N^-l 

X n / det[gf • xiwki)]--^ 

0<At<y<3 


where, following the same steps as taken in Ref. 
indices k and I in the set {!,... ,n}. 


one has now for each pair of different 


Vipki) 





(Dl) 


In the expression above, as a straight forward generalization of a 4-point fermionic Green’s 
fnnction, one dehnes su = {pk + and bki = yk± — yi± is taken in the k — I center of 
mass system. From here, following again the calcnlation of Ref.jd], one arrives for a given 
monomial of the expanded Vandermonde determinant at a strnctnre of, 

sm[al{OT)i\ sin[a;^(C>T)i] 


n 


nil 


l<k<l<n 


xGll 


kl 


2=1 




k(^i 


z{-ir] 
' EliPh 


r*+oo 


b + CxD 


32M 


-i \g(p{bki)J m 


dal 


3-2gi 1 . 

4 ; 2’ 


1 

2 ’ 


2 ’ 



ai 


a. 


where \pk\ = \pi\ = Pku and Ek = Ei = Eki, in each of the 2 to 2 -‘k — V— center of mass 
system. A 2n-point fermionic Green’s fnnction therefore is given as the prodnct of 2 
by 2-bnilding blocks that, np to their own kinematical variables, are identical to that of the 
starting expression o (as shonld appear clear in Appendix F, in order to be able to write 
the above generalization of ([T]), the resnlts of Appendix F have to be anticipated). 
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Appendix E: Integration of expression ([T]) 


The same steps as followed in ^ can be taken here to bring the second and third lines 
of ([I]) into the form, 


±( 


IQtiw? 

^2 


N 


N .0 

*=1 >^0 


da\ sin [k a.\{OT)i 


X / da2 sin [k 02 (C>T)i] (^34 I {a\a\Y 


0 , 

1 

2 ’ 


l-2qi 1 
4 ’ 2 


1 

2 ’ 


0 , 0 


For the sine function’s 

^/m2 Ep 


where k is here introduced as a shorthand for the constant 
arguments, no summation is understood over i, which is a distinguished index, and as usual. 


the overa 
of Ref. 


1 brackets denote the OAr(M)-average. Also, with respect to ([T]), formula 5.3.1(8) 


28j has been used. 


In the end, parameters are such, in the sine functions as well as in the Meiier special 
functions, the array of numbers { ar{i) } and { hg }, that formula 20.5.(7) of Ref. [29(] can be 
used twice so as to take the above expression into the form of. 


±( 


X 


■ 


N 


N 


nil-4-1) 


Qi] 


i=l 


/ y/32iNc Ep 




[jon 

9^{b) 


1-2 


^52 

■ CT4y 


I —14(0^1- 

Ep ')' 


1 , 1 , 


1 

2 ’ 

111 2qi+3 

2’ 2’ 2’ 4 ’ 


1 

2 

1 , 1 , 



can 


where the inversion formula 5.3.1(9) of Ref. 28|] has been used in order to express the result 
in a form suited to a |z| < 1 expansion of the Meijer function. 

Again, in the expression above, the array of numbers is such that 5.3.1(7) of Ref. 
be used twice so as to eventually take it to the much simpler form of, 

N 

nil-4-1) 


±( 


Ibvr^m^, 


N 


^2 


Wil 


i=l 


X 


V m2 y gip{b) 


1-2 


^03 


g^p{b) 

\/hl2iN^Ep^ 


\OT). 


1 3 -|- 2qi 
2’ 4 


(El) 


Note that the integration formulae of Refs. 


28 


29j have been extended to matrix-valued 


arguments. That this makes sense results from an available expression of Meijer’s special 
functions, suited to the case of \zi\ < 1. In the two previous forms, the Meijer function 
argument, Zi, is proportional to a quantity denoted by A in the main text, and defined in 
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(jH). Even at sufficiently large values of the coupling, = 10 — 20, for instance, this number 
is very small in the range of moderate sub-energies s. In that single quark model, taking 
a light dynamical quark mass of m ~ 5MeV, and a moderate sub-energy \/5 ~ 100 MeV, 
even at a (very) large non-perturbative mass scale of about 300 MeV for /i (corresponding 
to a very large estimate of the chiral condensate, for example), one has 


0.013 < A < 0.13, 


(E2) 


while the matrix elements of [(C^T)*]^ are themselves of order unity: For all i,j, in effect, 
one has < 1, and so are also all of the matrix element of the Tj, for all j = 1,..., iV; 

and thus, a fortiori, their product. 


For \z\ < 1, a Meijer function such as can be written as |28|, 




h=l 


bi, ^2, 63 j — ~ X 0-^2 [1 V bh — hi,... . .1 + hh — —d]E3) 


where 0^2 is a generalized hypergeometric series in z, with 


28|, 


0^2(1 bh{qi) - bi{qi), ..., =t=,..., 1 -F bh{qi) - b 3 {qi); -z] 


1 + Oih z + 0{z'^), 


(E4) 


where in flE3p . the prime on top of the product indicates that the pole value bj = b^ is 
omitted, and likewise the asterisk in the 0 F 2 generalized hypergeometric series (IE4I1 . indicates 
that the value bj = bh is omitted. This shows that for \zi\ < 1, the last expression above 
admits the expansion. 


±( 




N 


N 




h J^h 2 


1 + Oih Zi -\- Oi^zf )) 


i=l 


h=l 


Oiv(®) 


where (1E3I1 and (lE4p have been used, and where the variable Zi is the full (matrix-valued) 
argument of the Meijer function dehned in (01). The numbers and Oih are those given 
in dH]) and ([7]), respectively. Averaging over Ojv(M), the odd powers of {OT)i vanish, but 
there remain contributions of order z^, ^Jzi, Zi and Zi^JYi, which are on the order of Isxs, 
{OT)f, {OT)j and iOT)^ respectively, and are leading in view of the smallness of A. 

This illustrates that any order n in an expansion of the Meijer function is well dehned as 
an expansion in integer powers of matrices [{OE)i\, and this justihes the second and third 
expressions of this Appendix, or of ([3]) in the main text. 
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Appendix F: ‘OAr(M)-averaged Baker-Campbell-Hausdorf’ 

Upstream of ([T]), and leading to the second line of ([T]), one has an expression of jd], 


^ p-\-oo 

'/iCN^l[ dji e-4«p[(5i_5^) /jK.). (FI) 

i=l J-oo 


where fi{^i) = (exp [—and 6* as in (lAlOll . Then {N = D x {N^ — 1)), some 
convenient change of integration variables snch as dehned after flAlOll allows one to rewrite 
the hrst line of (1FT]1 as 14], 


2' V2 


/ d^ai / / d^Ui / d^U^e'^W, 


where, as thronghont the cnrrent paper, the brackets are here to mean an (9Ar(M)-averaged 
qnantity. The two exponentials of argnments aj-OT, (J = 1,2) are those leading to the two 
sine fnnctions that appear in the second line of ([T]), once integrations over the Uj, J = 1, 2 
are performed (level of expression flAlljl h 


Now, the point is the following. Whereas the 2nd line of flFlD displays a strnctnre of 
a hnite snm of hnite prodncts comprising N terms each, the 1st line doesn’t. This is 
of course due to the fact that both exponents are linear combinations of non-commuting 
SUc{3) generators (Cf. Baker-Campbell-Hausdorf formulae), so that in general, one has, 

N 

^iaj-oT ^ , J=l,2. (F2) 

i=l 

Now, it is a most important fact that under O^r(M)- integration, equality is restored in flF2F 
We will proceed by inspection of the hrst few non-trivial orders and check, 

N 

F”'‘°T,e = (n®”"‘°^’‘he . (F3) 

i=l 

At hrst non-trivial order (and besides the fact that odd powers of O trivially vanish under 
Ojv(M)-average) this is just, 

N N 

(1 + i ^ ai or,), e = (11 (1 + “) On)>,,e (F4) 

i=l i=l 
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At next order, one has to compare (simplifying to the case of iV = 2, so as to avoid too 
lengthy expressions), 

+ o'ja] (OT,OT, + OT^OTi) + 

to, 

f + 2 aWj {OT1OT2) + „ , 

that is, to compare the crossed terms. Within the notations introduced for Equation (ITOD . 
one hnds a result of i'^ajUj >e for both crossed terms, and equality holds. 

At order 3, one has to compare longer expressions. These match if and only if, 

((OTUOTUOTh + {OT^foTi)^^ = ((or),(or)2(or)i + = o. 

(F5) 

Now, besides the fact that, again, odd powers of O trivially vanish under OAr(M)-average, it 
is easy to check that the sole e-average guarantees flFSp . 

Fourth order is less trivial. Even at = 2, the corresponding expressions are quite 
cumbersome. At this order, the comparison of both sides of flF3p is between the 2 members, 

^{{oriOTl + {OT 1 OT 2 Y + OTiOTjOTi + (1 o 2 )))^q 

and. 

All other terms, comprising odd powers of OT\ or (9T2, give zero under full OAr(M)-average. 
In the 1st member, the terms with subscripts 1 and 2 exchanged just amount to twice the 
contributions of the hrst 3 terms. Then, resorting to explicit calculations one hnds. 


(or;or^>^g = (oTiorloro^g = ((OTior,)^)^ ^ 


N 


*J=1 




(F6) 


where, in the second line the SUc{S) constants dijk are dehned with respect to the values of 
indices i,j,k understood modulo 8. The two hrst equalities of flF6l) guarantee the equality 
of both sides of flF3p at fourth order since one has now. 


2 X 3 X 


4! 


2^,1 ^,2 


a ja 


J^J\2 


2 ! 


(F7) 
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Throughout this calculation use has been made of relations such as, 


(a«(©)a“(©)>e = <5'‘((a«)2(©)>e, 


(F8) 


met already at the level of (ITU]) , and which are easily derived out of (jUj) by inserting, between 
rotators, the relevant number of identity resolutions. In = N >< *1- 

In the case of two exponentials, relevant to fIFll) . exactly the same proof goes through. 
Calculations are even easier: At = 2 for example, it is immediate to see that by the 
discrimination of (ctjOj)^ into at any given order, the number of terms to be 

considered is diminished, and one gets in the same way, 

N 

^ia2-OT)^^ = ( Yl (F9) 

i=l 

That is, over exponentials of non-commuting S't/c(3)-generators, Ojv(®)-integration realizes 
the same simplification as Time-Ordering does on non-commuting field operators, 

j' (FIO) 


Therefore expression ([T]) is legitimate, with the accompanying remark that on the right hand 
side, the index i is not to be summed upon in expressions a.\ {OT)i and 
These are the two exponentials that, one step beyond, being integrated over aj from —oo 
to -I-CX3, together with another factor of l/oj (resulting from integration on V} in (lAlip b 
yield the two sine functions appearing in (IT]) . 


Relation (IF9P extends in a straight forward manner to the case of 2n-point Green’s 
functions and guarantees their non-trivial and interesting overall structures, as finite sums 
of finite products of Meijer special functions. 

As will be explained elsewhere, that (]F9]) has to hold true is indeed related to the structure 
of the Wick theorem expansion for Green’s functions. This structure must be the same 
in either perturbative and non-perturbative cases, and (IF9p guarantees that equivalence. 
Alternatively, this equivalence sheds some interesting light on the deep meaning of effective 
locality. 
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